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We examine the role that the gap dependence of the pairing interaction plays in the gap equation 
for a weakly coupled uniform superfluid of three-flavor massless quarks near the transition temper- 
ature T c . We find that the feedback effects on Landau-damped transverse gluons mediating the 
pairing interaction alter the gap magnitude in a way dependent on the color structure of the gap. 
We estimate corrections by these effects to the parameters characterizing the fourth order terms in 
the Ginzburg-Landau free energy and ensure the stability of a color-flavor locked state near T c . 
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I. INTRODUCTION 

It has been noted since the seminal work in Refs. 0, 0] that the quark-quark interaction in the color antitriplet 
channel is attractive and drives a Cooper pairing instability in quark matter in the limit of high density in which 
^C) ' the Fermi energy of the quarks dominates over the one-gluon exchange interaction energy. Because of the relativistic 
7-H , nature of high density quark matter, the color magnetic (transverse) force plays a dominant role in Cooper pairing. 
This interaction is nonlocal in time just like the electron-phonon interaction in ordinary superconductors UI3> but in 
contrast is long-ranged in the absence of static screening. This long-range nature gives rise to a different dependence of 
the pairing gap on the QCD coupling constant g from the BCS result 5j. At zero temperature, up to leading order in 
5, the logarithm of the gap arises from the dynamically screened magnetic force that involves Landau damped virtual 
gluons in a normal medium. Corrections to the gluon self-energy (polarization function) by the gap do not affect the 
1 logarithm of the gap up to subleading order in g at zero temperature This is because the gluon self-energy is 
modified by the gap significantly only at gluon energies of < T c . The influence of such corrections on the gap equation 
• near the transition temperature T c , however, has yet to be examined. 

In this paper we investigate the structure of the gap equation near T c by including the polarization effects of the 

■ color superconducting medium on exchanged gluons. We then estimate corrections thereby induced to the parameters 
fi |j characterizing the fourth order terms in the Ginzburg-Landau free energy of a weakly coupled uniform superfluid of 

massless three-flavor quarks. These corrections, divided by the weak coupling value, are of order g, in contrast with 
' the case of a short-range pairing interaction in which the corrections generally contain a factor proportional to the 
ratio of the transition temperature to the Fermi energy. We find that the polarization corrections keep the color-flavor 
locked phase the most stable just below T c . Throughout this paper, we consider a system of three-flavor (uds) and 
. \ three-color (RGB) massless quarks at temperature T and baryon chemical potential [i, and use units % = c = 1. We 

■ assume that the Fermi momentum is common to all colors and flavors. 



II. GAP EQUATION 



In this section we address the question of how the gap equation relevant in the weak coupling regime, in which the 
pairing interaction is induced by one-gluon exchange, is modified by the polarization effects of the color superconduct- 
ing medium. For this purpose we first consider a J p = + pairing state that is u<i-isoscalar and i?G-color antitriplet, 
since this is one of the simplest states that belong to a color and flavor antisymmetric channel with J p — + . This 
channel has a common transition temperature in the limit of weak coupli 
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where e 7 /n = 0.5669 . . . and b = 2567r 4 (2/3<7 2 ) 5 / 2 . For the isoscalar pairing state of interest here, a nonzero excitation 
gap d(k) is open below T c for quasiparticlcs having either flavor u or d and either color R or G. Here k is the 
momentum associated with the relative coordinate of a quark Cooper pair. In the case in which only the modification 
of the one-gluon exchange force by a normal medium is included in the random-phase approximation (RPA) 0, |jj 
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and the normal state Hartree-Fock correction to the quark propagator is ignored, the gap equation reads jToL ITT| 

d(k) = yd 3 g[^ T (i?(q)-^(k),q-k)+^ T (i?(q)+ii;(k),q-k) 
+D L {E{ci) - E(k), q - k) + D L (E(ci) + E(k), q - k)] 



where 



is the excitation energy, 
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D T (p) ~ Re 
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with the Debye screening mass 
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characterize the transverse and longitudinal parts of the gluon propagator D(p) in the Landau gauge as 

DfAP) = -5 a p[P^D T { P ) + P^D L (p)], (7) 
with the transverse and longitudinal projection operators 

T r PiPj ^J 1 / \ 

Mi = "ij ~ I 12 ' M30 = °0« = °»0 = 0) (8) 



(9) 



Expressions Q and JSJ are approximate in the sense that they are available in the regime po <§; |p| <C /x/3, but duly 
allow for the Landau damping of transverse virtual gluons and the Debye screening of the longitudinal force in a way 
sufficient to describe the exact form of the logarithm of the gap magnitude up to subleading order in g |l0(. Note 
that the Landau damping provides an effective infrared cutoff in the transverse sector, ~ (nm 2 D \po\/4:) 1 / 3 , which in 
turn plays a dominant role in determining the pairing gap. 

At T = 0, the solution to the gap equation J5J is known as 



d(k) = -biMe-*& sin(^), 



where g = g/3y2w, and 
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The factor sin(^x) ensures that the gap is appreciable only for momenta k close to the Fermi surface. The exponential 
term and the sinusoidal x dependence arise from nearly static, Landau-damped magnetic gluons that mediate the 
long-range part of the magnetic interactions, while both the higher frequency magnetic gluons and Debye-screened 
electric gluons play a dominant role in determining the pre-exponential factor. 



3 



For later comparison with the case near T c , we write down the equation for the magnitude of the gap on the Fermi 
surface, dp = d(\k\ = /i/3), which can be derived from Eq. (J^Jl as [l(| 
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Here the cutoff 5, obeying dF "C S <C mp>, is chosen so that d(|k| > S) is vanishingly small, and b' — 2 ln(&/i/3<5). The 
term associated with hi 2 (28 /dp) comes from soft Landau-damped magnetic gluons, while the term associated with 
ln(2<5/dp) comes from nonstatic magnetic gluons and Debye-screened electric gluons. 

The overall coefficient of the n/g 5 in the pre-exponential factor in Eq. (|10|l is correct up to a factor of order unity 
since the quasiparticle wave function renormalization ignored here results in a factor exp[— (n 2 + 4)/8] |12|. which 
appears also in the weak coupling expression for T c , Eq. Q. (This renormalization affects the sinusoidal x dependence 
only through a factor of order g 2 x.) On the other hand, polarization by the color superconducting medium, which gives 
rise to the gap dependence of the pairing interaction, provides even higher order corrections to the zero-temperature 
gap @. 

Near T c , the momentum dependence of the gap can be set equal to that at T = 0, as in the usual BCS case |10|. 
Consequently, 
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d(q,T) = d F (X)sm(gy), 
26/x/3 



|q|-/i/3|+£(q,T = 0)_ 
Then, the gap magnitude dp can be determined by expanding the gap equation @ up to 0(d 3 ) as 
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/x/3| 



(15) 



where we have noted that the momentum region, k w q, contributes dominantly to the integral in Eq. J2Jl. We thus 
obtain 



dp = 1 - 



izg T — T c 
~2 t7~ 



16^T 2 F 



Fh 



(16) 



with the zeta function £(3) = 1.2020 . . .. Here the coefficients affixed to dp and dp include the leading contributions 
with respect to g, and we have noted that at T = T C , 
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is satisfied. [In the absence of the quasiparticle wave function renormalization, the solution to Eq. (|17f) reproduces 
expression (QJ except for a factor exp[— (n 2 + 4)/8].] The solution to Eq. (fTB)l reads 



dp = 



8tt 2 T 2 T - T r 



-,1/2 



7C(3) T c 



(18) 



As it should, Eq. (|16|l is the same as the known result obtained from the Ginzburg-Landau theory (see the next 
section). 

We now introduce the effect of the color superconducting medium on the gluon propagator within the RPA. The 
normal gluon propagator characterized by Eqs. J2J and (J5J is modified by the pairing gap, in a way dependent on a, 
as 131 



D%(p) ~ Re 
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IPl 2 + (m a M ) 2 f(p) - iirmlp d(V^m D /2 - |p|)/4|p| 



(19) 
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and 



D£(p) ~ IpP+m^-atm^Mp)- (20) 
Here only leading corrections by the pairing gap have been retained, and 

C 0, a = 1,2,3, 

(™li) 2 = { 9 2 K T d 2 F , a = 4,5,6 ,7, (21) 

[ {4/3)g 2 K T d 2 F , a = 8, 

are the Meissner screening masses |14| with the stiffness parameter in the weak coupling limit: 
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where 
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is the ideal gas density of states at the Fermi surface. f(p) and h(p) are the dimensionless positive definite functions 
that characterize the 0{d 2 ) corrections to the transverse and longitudinal parts of the irreducible particle- hole bubble 
[Tflj . These functions reduce to unity in the limit of p — > 0, while decreasing to zero with increasing |p|. We note 
that the Landau damping term, corresponding to the energy-dependent term in Eq. I|19(l . does undergo corrections 
by a factor of 1 + 0(cP/|p| 2 ), but they lead to higher order corrections to the gap equation in T c /p as compared with 
those coming from the Meissner term (m^) 2 / (p) . 

Expressions H19|) and H2U|) are the straightforward extension of the normal medium forms (@J and JSJ to the case 
of the color superconducting medium. These expressions for Df and retain consistency with the transverse and 
longitudinal sum rules obeyed by the static, long-wavelength gluon propagator in normal quark matter p^j | and in 
color superconducting quark matter [l4|. 

It is important to note that the leading feedback effect of the pairing gap lies in the magnetic sector. This is partly 
because the corrections to the gluon propagator are roughly of order m 2 M / '(m^lpol) 2 ^ 3 in the magnetic sector, while 
of order m 2 M /m 2 D in the electric sector, and partly because the gluon energy range \p<y\ <C mo plays a dominant role 
in the gap equation. In deriving such corrections in the magnetic sector, we first expand Dt, Eq. Ijl9|l . with respect 
to m^, and then substitute into /(p) the form relevant near T c : 

6 °° f 1 1 — x 2 

/(P) = 7C(3) §/ dX (s + 1/2)[4( S + 1/2)2 + Q P \ X /2*T D )*] " (24) 

This form is identical to that encountered in the usual BCS case 0| since in both cases a pairing gap is open for 
quasiparticle momenta so close to the Fermi surface that quasiparticles and quasiholes having momenta ki and k2 
with + k 2 = p and |ki| ~ |k 2 | ~ n/i dominate the 0{d 2 ) corrections to the particle-hole bubble. In the London 
limit (|p| — > 0), one can set / — 1, while in the Pippard limit (|p| — > oo) / behaves as oc l/|p|. As we shall see, the 
Pippard regime (|p| > 2irT c ) is as important to the gap equation as the London regime (|p| < 2irT c ). 
The gap equation modified by the leading feedback effect can be obtained by replacing by 



D T {p) ~ Rc 
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|p| 2 - iirm 2 D p 6(y/^m D /2 - |p|)/4|p| 



(m^) 2 /(p)Rc 



1 



|p| 2 - mm 2 D p e{^m D /2 - |p|)/4|p|] 2 



(25) 

in Eq. J2J. The modification associated with m%j works only for a = 8. This is partly because the Meissner masses 
vanish for a = 1-3 [see Eq. (|21f> ] and part ly because the contributions of gluons of a — 4-7 to the gap equation vanish 
due to the color structure of the gap |10|. Near T c , the gap equation at |k| = /Lt/3 thus reads 

d F = I1-— ' " c \d F - '^ K ~ll d% 
! 2 T c J 16ttT 2 f 

2 r & 

+0(d%), (26) 
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where 



F(q) - J 



d(cos 8) — 



1 — cos 8 



) 3 + {irm 2 D 



|q| - M /3|/8\/2|q|3)20(0F mD /2 - y/2\q\ y/1 - cos ( 



/(!/) + !/ 



4f(w) 



_ (27) 

with y = |q|\/l — cos8/^/2irT c . Here, is the angle between k and q, we have replaced |k| by /i/3 in the gluon 
propagator, and the factor 1/12 comes from the color vertex part of a = 8. Note that the leading term by the 
feedback effects in Eq. I|26[) is of third order in cLf and hence does not affect T c . 

We proceed to examine the leading corrections to the gap equation due to the color superconducting medium. For 
this purpose, it is useful to divide the gluon momentum |p| ~ v2|q|\/l — cos6> and energy po ~ ||q| — /x/3| into several 
regimes shown in Fig. 1. The boundary y = 1 corresponds to the gluon momentum |p| = 2irT c . For y < 1 (London 
regime) and y > 1 (Pippard regime), f(y) has the following forms: 



}^ a n y 2n , a n 

n=0 



24 (-l)"C(2n + 3) 
7C(3) (2n+l)(2n + 3) 



and 



E 



l 



2(s + 1/2) 2 y 2 



tan 



Y _ 2-(2n+3) 
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2(s + l/2) (s + l/2)2/ 



(28) 



(29) 



At y = 1, as can be seen from Eq. 127|) , the effective infrared cutoff due to the Landau damping of transverse gluons is 
nonnegligible when po is larger than 327r 2 T c 3 /m| ) . Another important scales are pq — 7rT c , above which tanh(po/2T c ) 
in Eq. i|2ti[) approaches unity exponentially, and |p| = \/Ttmz)/2, above which transverse gluons no longer undergo 
Landau damping. We note that the gap can be regarded as flat in regions a) and c). We also note that the timelike 
regime (po > |p|), in which the gluon propagator is not described well by the form (|25|) . can be safely ignored since in 
this regime no feedback corrections to the gap equation occur up to leading order in Tjji. In fact, the gap corrections 
to the gluon propagator vanish like ~ d 2 jp\ with increasing p§, as in the T — case 
At 2irT c < |p|, we can use expression (12911 for f(y). From this expression, we obtain 



f(y) + y 



12 



df(y) 

dy 7C(3)y ; 



;Hy), 



(30) 



with 



Hy) = E 



2 _ x y 1 
y 2(s + 1/2) + s + 1/2 



The large y asymptotic behavior of f(y) + ydf(y)/dy is then 



f(y) + y 



12 



df(y) , 

dy ~ 7C(3)2/ : 



Iny. 



(31) 



(32) 



This behavior is different from cx which is followed by f(y). As a result, transverse gluons of momenta near T c 
rather than near the Pippard limit are essential to calculations of the feedback effect in the magnetic sector. This 
is a contrast to the case of the weak coupling limit in which transverse gluons of momenta large compared with T c 
dominate the pairing interaction since for such momenta, the factor jti^po/IpI characterizing the Landau damping in 
the propagator (JIJ is sufficiently small that the pairing interaction remain essentially long ranged. 

Using Eqs. I|27|l and l|29|) . we calculate the contribution to the gap equation l|26|l from regions c) and d). The 
momentum range covering these regions corresponds to the range of I 
leading order in g, the result from region c) reads 



satisfying -1 < cosfl < 1 - 27r 2 T c 2 /|q| • Up to 



7T[3^ 3 - 28C(3)/(l)]5 2 (m^ 8 ) 2 d 



112C(3)m 2 1 

while the contribution from region d) is of higher order in T c j \i and thus can be ignored. 

At |p| < 2ttT c , where expression (12811 is available for f(y), f(y) and hence f(y) + ydf(y)/dy are almost flat. 



(33) 



In 



-F(q), this momentum range corresponds to the range of 8 satisfying 1 
the gap equation 12(j|) from region a) becomes, to leading order in g 1 



■Kg* 



4m 2 D 
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71=0 



27r 2 T c 2 /|q| 2 < cosfl < 1. The contribution to 



(34) 
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FIG. 1: The energy-momentum regimes of exchanged magnetic gluons. 



while that from region b) is of higher order in T c / fj,. It is remarkable that the term l|34|) is comparable to the term 
(|33[l . This suggests that transverse gluons of momenta below and above 2irT c are equally important to the feedback 
effect. We also note that these terms are of order g 2 dp/T 2 and thus suppressed only by one power g with respect 
to the term proportional dp in the gap equation (|16|l in the weak coupling limit. This is a contrast to the case of 
a short-range pairing force in which the leading correction to the third order term due to the superfluid medium is 
suppressed by one power T c //j, [16J. We remark that the scale of the gluon energy dominant in the gap equation is of 
order 327i- 2 T c 3 /m|j. Since this is much smaller than the typical momentum scale ~ 2ttT c , we can safely take the static 
limit of the modification due to the color superconducting medium. 

We finally rewrite the gap equation 1|26[) in such a way that the feedback term has a coefficient up to leading order 
in g, i.e., by combining the contributions (|33[1 and 1|34|) from regions c) and a). The result is 

*-('-?^)*-&' + Z ^* + °<*>- 

where C = ?r 3 /63C(3) = 0.409434 .... We thus find that the leading feedback effect near T c acts to increase the gap 
squared of the isoscalar pairing state by a factor of (1 — Cg/2)~ 1 . This is due to the fact that the feedback effect 
manifests itself as Meissner screening of the color magnetic force of color index a = 8; this force is repulsive in contrast 
to the attractive case of a = 1-3 dominating the pairing interaction. We remark that in the T = case in which 
expansion of the gap equation with respect to the gap magnitude is not valid, the terms associated with the logarithm 
of the gap in the gap equation l|12|) mainly determine the gap magnitude. In this case, the feedback effect provides 
corrections beyond these logarithmic terms [||. 

For the purpose of calculating corrections to the Ginzburg-Landau parameters in the next section, it is instructive to 
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repeat the above calculations for a color- flavor locked (CFL) state, one of the J p = + , color and flavor antisymmetric 
pairing states. In the CFL state, all three flavors and colors are equally gapped in such a way that the pairing gap 
between a quark of color a and flavor i and a quark of color b and flavor j is characterized by K(5 a iSbj — S a j Su). For 
the on-shell gap on the Fermi surface, kf, the gap equation near T c can be written in the form similar to Eq. (|35[l as 

ngT-T c \_ 7C(3)^ 3 21C(3)Q 2 „ 3 , ^. 5 



kf = { 1 - f -fr) ^ - - -g^** + w). (36) 

Here we have used the Meissner masses in the CFL state [l4|, i- e -i ( to m) 2 = ^9 2 Ktk f for a = 1-8. In the limit of 
m M ~ * 0' Eq. l|3*6*)) is equivalent to Eq. (104) in Ref. We thus find that the leading feedback effect acts to reduce 
the gap squared by a factor of (1 + 3Cg) _1 . This reduction stems from the fact that Meissner screening of the color 
magnetic force takes effect equally for a = 1-8. Note a contrast with the case of the isoscalar pairing state in which 
the feedback effect acts to increase the gap magnitude. 

III. GINZBURG-LANDAU FREE ENERGY 

We proceed to derive the Ginzburg-Landau free energy of a weakly coupled uniform superfluid of massless three- 
flavor quarks from the gap equation near T c as examined in the previous section. Instead of focusing on the isoscalar 
and CFL pairing states, it is convenient to construct the Ginzburg-Landau free energy for a more general color and 
flavor antisymmetric channel with J p = + as in Ref. [TlJ. This is because all states belonging to this channel has a 
common value of T c , which reduces to Eq. Q in the weak coupling limit. This channel is characterized by a complex 
3x3 gap matrix, (d a )j, in color-flavor space where a (i) is the color (flavor) other than two colors (flavors) 
involved in Cooper pairing. This gap is defined on the mass shell of the quark quasiparticle of momenta on the Fermi 
surface, and thus reduces to 5 a B5i S dp in the isoscalar state and to 6 a iKF in the CFL state. For (d )j, one can write 
down the thermodynamic potential density difference Af2 = Q s — VL n between the superfluid and normal phases near 
T c as [ll| 

a a ab 

Here each term is invariant with respect to U(l) global gauge transformations and color and flavor rotations. 

In evaluating the coefficients in Eq. (|37() by including the leading feedback effect, it is convenient to integrate the 
gap equations l|35|l and l|36|l and then map the results onto Eq. (|37|) . We thus obtain 

a = 4N{n/3)hx(^r\ (38) 

A ^H(i + y^)%/3). (39) 
A = 8ra i ( 1 -T C1 ')^ (40) 

which reproduce the known relation (3\ — /3 2 in the weak coupling limit [Tlj . We find that up to leading order in 
g, the polarization effects of the color superconducting medium give rise to only 0(g) corrections of the coefficients 
Pi and /?2- It is nonetheless important to note that these corrections to Pi and P2 work in the opposite directions 
in such a way as to decrease and increase the gap magnitude in the CFL and isoscalar states, respectively. We can 
thus conclude that whether or not the leading feedback effect acts to reduce the gap magnitude depends on the color 
structure of the gap. 

IV. PHASE DIAGRAM 

We turn to the construction of the phase diagram in the space of the parameters characterizing the fourth order 
terms in the Ginzburg-Landau free energy derived in the previous section. This phase diagram can be obtained by 
minimizing the thermodynamic potential difference AJ7 with respect to (d a )j for various values of Pi and P2 The 
result, exhibited in Fig. 2, is the same as Fig. 1 in Ref. 01 except that the present values of pi and P2 given by Eqs. 
(13911 and (|40|l include the O(g) feedback corrections. This figure shows that even for such values of Pi and P2, as long 
as g is sufficiently smaller than unity, the CFL phase, which generally satisfies 

d| J -d G = d^-d B = d| ; -d i? = 0, |d fl | 2 = |d G | 2 = |d s | 2 , (41) 
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FIG. 2: Phase diagram near T c , exhibiting regions where the 2SC and CFL phases are reached by a second order transition 
as well as where superfluid states are reached by a first order transition since the overall fourth order term in the Ginzburg- 
Landau free energy (13 7H can be negative. The parameters f3\ and fa are the fourth order coefficients in the Ginzburg-Landau 
free energy. The cross denotes the weak coupling limit, and the circle denotes the result including the polarization effects of 
the color superconducting medium with g = O.f . 

is still favored over the two-flavor color superconducting (2SC) state fulfilling 

d R || d G || d B - (42) 
Note that the 2SC state contains the isoscalar state analyzed in SeclTTl 

V. CONCLUSIONS 

We have examined the role played by the gap dependence of the pairing interaction in the gap equation for a 
weakly coupled uniform superfluid of three-flavor massless quarks near the transition temperature T c . The corrections 
induced by this role to the parameters characterizing the fourth order terms in the Ginzburg-Landau free energy result 
in an increase of j3\ by a factor of 1 + HCg/2 and an decrease of /?2 by a factor of 1 — 15Cg/2. The magnitude of these 
changes comes from the Meissner screening of the color magnetic force that dominates the interaction between quarks 
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in the weak coupling regime as an essentially long-range force. We thus see a contrast with the case of a short-range 
pairing interaction in which the corrections are suppressed by one power T c //x. We also note that the changes in /3± 
and @2 are in the direction of increasing the gap magnitude of the 2SC state and decreasing that of the CFL state. 
This direction reflects the fact that not only the Meissner screening of the color magnetic force but also the color 
indices of magnetic gluons dominating the pairing interaction depend on the color structure of the pairing gap. We 
have finally found that the feedback corrections, as long as g -C 1, keep the color-flavor locked phase the most stable 
just below T c . 

The present result for the parameters /3i and P2, mainly through its effect on the gap magnitude near T c , provides 
a way of studying strong coupling modifications on the previous weak coupling calculations based on the Ginzburg- 
Landau theory with fli = f3 2 - Those calculations include the phase diagram iTll as discussed in Sec. Eland its 
extension to nonzero quark masses [T^| . responses to rotation and magnetic fields [La. 1 19Ll20| . and fluctuation- induced 
first order transition |2l|. Qualitatively, however, no significant changes are expected. Nonetheless, the tendency that 
Pi > Pi-, if remaining at low densities, could be significant for the normal-super interfacial energy and the interaction 
between widely separated magnetic vortices for CFL quark matter. This is because both quantities, which are sensitive 
to the ratio between j3i and 02 [Til 12(1 ] , control the criterion of whether or not the CFL state can allow magnetic 
vortices to form. 
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